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Abstract 

We quantize the elastic modes in a plate. For this, we find a complete, 
orthogonal set of eigenfunctions of the elastic equations and we normalize 
them. These are the phonon modes in the plate and their specific forms 
and dispersion relations are manifested in low temperature experiments 
in ultra-thin membranes. 

1 Introduction 

Nowadays, the devices used in many high-sensitivity applications reached such a 
level of miniaturization, that the wavelength of the quantum quasiparticles used 
in their modeling is comparable to the dimensions of the device. The examples 
we are most familiar with are the ultra-sensitive electromagnetic radiation de- 
tectors. In a very general way of speaking, such detectors consist of some thin 
metallic films, a few tens of nanometers in thickness, deposited on a dielectric 
membrane. The dielectric membrane is usually made of SiN 2 , and has a thick- 
ness of the order of 100 run. To reach the level of sensitivity and speed required 
for applications, these detectors have to work at sub-Kelvin temperatures and 
at such temperatures the dominant phonon wavelength is comparable to the 

devices' thickness [TJIlElHISlEldliE]. 

To describe the thermal properties of such membranes or detectors, the 
elcctron-phonon interaction, or in general any interaction of phonons with im- 
purities or disorder in the membrane, we have to know the phonon modes in 
the membrane. For this, we have to find the eigenmodes of the elastic equations 
and quantize the elastic field. 

For infinite half-spaces, the quantization has been carried out by Ezawa |10j . 
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Figure 1: Typical plate, or membrane, like the ones used to support mcsoscopic 
detectors. The thickness d is of the order of 100 nm and I 3> d. The plate 
surfaces are parallel to the xy plane and cut the z axis at ±d/2. 



1.1 Elastic eigenmodes in plates with parallel surfaces 

The elastic eigenmodes in plates with parallel surfaces have been studied for a 
long time, mostly in connection with sound propagation and earthquakes (see for 
example Ref. [H]). To introduce these modes, let us consider a plate of thickness 
d and area I 2 , with / ^> d. The two surfaces of the plate, or membrane, are 
parallel to the (xy) plane and cut the z axis at ±d/2 (see Fig. [TJ. Throughout 
the paper we shall use V for the volume of the membrane (or in general of 
the solid that we describe-see Section 12. lj) and dV for its surface. We shall 
assume that I is much bigger than any wavelength of the elastic perturbations 
considered here. The displacement field at position r is going to be denoted by 
u(r) or v(r). The unit vectors along the coordinate axes are denoted by x, y, 
and z. 

The displacement fields obey the dymanic equation 

pd 2 u l /dt 2 = c ljM djd k ui , Vi = 1, 2, 3. (1) 

Here, as everywhere in this paper we shall assume summation over repeated 
indices. Assuming that the medium is isotropic (see for example Ref. [11] for 
the constraints on the tensor [c]), Eq. {T]) is reduced to d 2 Ui/dt 2 = djPij, with 
Pij defined as in [10], p i3 = p~ 1 c lj kidkui = (cf - 2cf )(dku k )8 ij + c 2 (d l u j + djUi), 
and q, ct the longitudinal and transversal sound velocities, respectively. 

Introducing the operator L (we shall use tilde for operators and hat for unit 
vectors) by Lu = p~ 1 CijkidjdkUi — c^graddivu — c^curlcurlu = cfV ■ V • u — 
C( V x V x u, the wave equation ([1]) becomes 

d 2 u/8t 2 = Lu. (2) 

The surface is free, so the stress should be zero there. This amounts to the 
boundary conditions [TTJ [TU] 

PijUj = 0, on dV, Vi = 1, 2, 3, (3) 

where n is the unit vector normal to the surface, of components n\,n2, and 
n3~we shall use this notation throughout the paper. 
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Applying Eqs. ^ and © to the plate we obtain the elastic eigenmodes, 
which are classified in three groups, according to their symmetry or polarization 
direction: the horizontal shear (/i), the symmetric (s) and the antisymmetric 
(a) waves. All these waves are propagating (or decaying, if the wave-vector 
is complex) along the membrane and have a stationary form in the direction 
perpendicular to the surfaces. The h wave is polarized parallel to the surfaces 
and perpendicular to the propagation direction. The s and a waves are su- 
perpositions of longitudinal and transversal waves, polarized in a plane that is 
perpendicular to the surfaces and contains the propagation direction. The dif- 
ference between the s and the a waves comes from the fact that the displacement 
field along the z direction is symmetric for the s wave and antisymmetric for 
the a wave, while the displacement field along the propagation direction is anti- 
symmetric for the s wave and symmetric for the a wave (see below). Explicitely, 
the three types of modes are: 



u,j = (k|| x z) • N h cos 



•mix 



z — 



,<(k u .r- W *) j m = 0,1,2,. 



(4a) 



N s \ z ■ fen 



+ky ■ ikt 

Xe t(k||T-U,t) i 

N a \ z ■ fc„ 



+k|| • ik t 



>*(k|| -r-urt) 



-2k t k t cos (kt^j sin (kiz) + [kf - fcjj] cos (kl^\ sin ( fc * z ) 

2fc^ cos (kt-^J cos (kiz) + [kf — fcjj] cos (kiT^J cos (fa z ) 
t) ) 

^)cos(fc/z) - [A: t 2 - fey] sin(fc;^) cos(fc t z) 



(4b) 



2k t ki sin(fc t -) 



2fc| sin(A: t |)sin(fc/z) + [k\ - kf\ sin(fc^) sin(k t z) 



xe 



(4c) 



where k|| is the unit vector along the propagation direction and k|| = k» ■ kii . In 
the s and a modes, the wavevector in the z direction takes two values, kt and 
ki, one corresponding to the transversal component of the mode, the other one 
corresponding to the longitudinal component. The constants Nh, N s , and N a 
are the normalization constants, which will be calculated in Section [3] 

The components of the wavevectors, kt, ki, and fc|| obey the transcendental 
equations [llj . 



and 



tan(|fc t 
tan(|fc; 



tan(|fc t ) 



Akik t kf, 



r 2]2 



[k 2 t-kW 



tan(|fc ; ) 

for symmetric and antisymmetric waves, respectively. 



(5a) 



(5b) 
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All the solutions of the elastic equation © are given by equations (f4"aj) . (|4b|) . 
and (|4cp . with to taking all the natural values, 0, 1, . . ., whereas k t and ki are the 
solutions of Eqs. (|5a|l and (|5b>|) . fen can be a complex number, but the complete, 
orthogonal set of phonon modes that we shall use in the quantization of the 
elastic field are the ones with running from to oo 

The paper is organized as follows: in section l2~Tl we show that L is self-adjoint 
even when applied to the displacement field of an elastic body of arbitrary shape. 
Therefore, we can form a complete, orthonormal set of its eigenfunctions. 

The fact that the elastic modes ((4]) corresponding to different quantum num- 
bers are orthogonal to each other is proved in Section 12.21 based on the her- 
miticity of L. 

The normalization constants are calculated in Section [3] and the formal pro- 
cedure of quantising the elastic field is presented in Section SJ 

We apply this formalism elsewhere to calculate the phonon scattering in 
amorphous, thin membranes. 

2 Orthogonality and completeness of the set of 
elastic eigenmodes 

2.1 Self-adjointness of the operator L 

We shall prove the self-adjointness of L on an arbitrary volume V. We assume 
that V has the smooth border, dV. The operator L acts on the Hilbert space 
TL which consists of the vector functions defined on V, which are integrable in 
modulus square. The scalar product on TL is defined as usual, 

(v|u) = f v^r) -u(r)d 3 r, (6) 
Jv 

and the norm is ||u|| = (uju) 1 / 2 . The domain of L, denoted by T>(L), is formed 
of such functions u(r) £ TL, so that iu(r) exists and is contained in TL. More- 
over, the functions from T>{L) should obey the boundary conditions ([3]). Ezawa 
showed that L is hermitian [TU]. We will show now that L is also self-adjoint. 

First of all, for a formal treatment, in this section we will understand the 
derivatives in a generalized sense. If f(r) is an arbitrary function on V, then 
dif is defined as the function that satisfies 

[ g(r)d i f(r)d 3 r= [ !{v)g{v) ni d 2 r - [ f(r)d i9 (r) d 3 r , (7) 

JV JdV JV 

for any function g(r) of class C 1 (V^) (i.e. g(r) is derivable, with continous first 
derivatives on V) and all the integrals on the right-hand side of Eq. ([7]) exist 
and are finite. 

Returning to our operator, let us first note that the space T>{L) includes the 
space of functions twice derivable, with continous, integrable, second derivatives, 
C 2 (V), which is dense in TL. Therefore T>(L) is also dense in TL, and L is a 



4 



symmetric operator. Then wc define the adjoint operator U and its domain 
V(L'). For this, let v be a function in TL, so that (v|Lu), is a continous linear 
functional in u G D(L), i.e. there exists an M v > for which 

(v|Lu) <M v ||u|| (8) 

for any u G T>(L). By the Riesz-Frechet theorem [T^], there exists a v* G TL so 
that (v|Lu) = (v*|u) for any u G T)(V). The adjoint operator is defined by 
the relation ZJv = v*, for all v that satisfy ©. The functions that satisfy © 
form the domain of L>, denoted T>{L^). Since L is hermitian, if v, u G T>(L), then 
(v|Lu) = (Lv|u) is a linear functional, so any function from T)(L) is included in 
V{D). Therefore we can write V{L) C V{D). To prove that L is self-adjoint, 
we have to show that also T>{U) c T>(L), so in the end T>(JJ) = T>(L) and 

£ = £t. 

For this, let us take v G T>(U) and u G T)(L). Integrating by parts we get 

p(v|Lu) = / v*c tjk idjd k uid 3 r (9) 
Jv 

{d j v*)c i:jk in k u l d 2 Y+ I (d k djV*)cijklUid 3 r 

dV JV 

{djV*)c l3k in k uid 2 Y + p(Lv\u) = p{Dv\u) , (10) 

dV 

where we used Cij k irijd k ui = on dV (Eq. [3]), and the simplified notation 

p(Lv\u) = / {ci jk idjd k v*)uid 3 r = / (ci kj id k djV*)uid 3 r , (11) 
Jv Jv 

although v is not necesarily a function in T>(L). The last equality in Eq. 

is obtained using Cij k i = Cji k i = Cijik = c k uj [H] and permuting the partial 

derivatives. But Eq. (jTOJ) means that 

({L - L f )v|u) = p- 1 [ d jV *n k c l0k iuid 2 v (12) 

JdV 

for any u G T>{L). If (djV*)n k Cij k i is not identically zero on 9V\ then we can 
find u so that the surface integral in (fT2"l) is different from zero. This implies 
that ((L — L')v\u) ^ 0, so v ^ on a set of measure larger than zero in the 
interior of V, denoted as V°. In such a case, we can find a nonempty compact 
set S C V° and a function, u', which is twice derivable, zero outside S, and 
satisfies 

<(Z-Zt)v|u')^0. (is) 

Since u' is zero outside S and S is a compact set in V°, this means that both, 
u' and any of its derivatives are zero on dV; therefore u' G T>(L). Now, u'(r G 
dV) = 0, implies that 

djV*n k c ijk i(ui + u' l )d 2 r = / djV*n k c ijk iuid 2 r (14) 

dV JdV 
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whereas (fT5)l implies that 

((X-£t)v|(u + u')>^((£-Xt)v|u>. (15) 

Since Eq. fT2]) should be valid for any function in T>(L), including u and u + u' 
and Eq. p4[) is true by construction, this implies that (fT5")) is a contradiction; 
therefore (djV*)nkCijki = on dV. Moreover, by the definition of T>(U) and 
the Riesz-Frechet theorem, cikjidkdjUi = pL^v G H. Therefore v S T>(L), so 
T>(U) = T>(L) and L = L>, i.e. the operator L is self-adjoint on an arbitrary 
volume, V. 

Going back to the plate with infinite lateral extension, we have to find from 
the wave functions of the form (0} a complete, orthonormal set. We can do that 
by using the operator k|| = i(d x + d y ), which is also self-adjoint when acting 
on wave-functions defined on a plate with infinite lateral extension, or on a 
finite, rectangular plate, with periodic boundary conditions at the edges. Since 
L and k|| commute and they are both self-adjoint operators, if we find all the 
cigenfunctions common to L and k||, then we have a complete set. But these 
functions are given by Eqs. (@|, with real k», and k t and ki satisfying ([5]). 



2.2 Orthogonality of the elastic eigenmodes 

Now we study the orthogonality of the elastic eigenmodes of the plate. For 
this, we write the the functions that appear in Eqs. (pia| - (|4"cl) in general as 
u k||,fc t ,o-(r) = Ufc t Cr (2;)e lk ii' r , where we separated the x and y dependence of the 
fields from the z dependence and we disregarded the time dependence. By a 
we denote the "polarization" h, s, or a. We shall not use ki explicitely in the 
notations below, since it is determined implicitely by k t , kn, and Eqs. ([5]). The 
operator k|| has eigenvectors of the form g(z)e lk H' r , where g(z) is an arbitrary 
function of z while the eigenvalue k||, perpendicular on z, has real components 
along the x and y directions. We analyse the common set of eigenvectors of L 
and k|| , so in what follows we shall consider only the cases with real wave- vectors 
kii . First we observe that 



,b/2 

(u k|l ,fc t , CT |u k ; k > „,) = 27T(5(k|| -kji) / u[ t Jz)u k> ,(z)dz, (16) 

J-b/2 



so we are left to study the orthogonality of functions with the same k||. For 
simplicity we choose k|| = x • fcy , so the h waves are polarized in the y direction 
and the s and a waves have displacement fields in the (xz) plane. Since the 
displacement field of the h waves are perpendicular to the displacement fields of 
the s and a waves, any h wave is orthogonal to any s or a wave. Similarly, for 
the same k|| , the displacement fields of any of the s and a waves, although in 
the same plane, are orthogonal to each-other due to their opposite symmetries. 
We conclude that 

,6/2 

/ U L A z ) u K,<T>( z ) dz = o, 

J-b/2 
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for any k t and k' t , if a ^ cr'. Therefore (uk||,fc t , CT |u k ' k' t ,a') #(k|| - k^)6 a ^. 

We are left to show that (uk||,fc t ,a-|uk||,fc£,a-) ^ <5fc t .fc; which follows simply 
from the hermiticity of the operator L. We calculate the matrix element 

(uk||,fc t , -|-£'Uk || ,fej, -) = -^£ |hfe ; !CT (u k „,fc t , ff |u k|hfc / !0 .) 

= (£u k||ifctjCr |u k||j fc< i0 .) = -Wk| h fc t ,o-( u k| 1 ,fc t , -|uk |h fe;, '} 

But since for given k|| and cr, the eigenstates of L are not degenerate, (u k||) fc tiCr |u k fc/ )Cr ) = 
0, unless k t = k' t and, of course, ki ~ k[. This completes the proof and 

(ukn,fe tl0 -|uk' k' t ,<r>) oc <W5(k|| - kj|)4 tifc / . (17) 

In the next section we calculate the normalization constants. 

3 Normalization of the elastic modes 

From this point on wc shall assume that the plate has finite lateral extensions of 
area A = I x I and the eigenfunctions obey periodic boundary conditions in the 
(xy) plane. The volume of the plate is V = Ad. As a consequence, the values 
of k|| are discrete and the scalar product (JT7J) should give 

(u k|| , fcti(T |u k / ifc /^) = S a , a >S kll .^5 kt . k ' t . (18) 

For horizontal shear waves, the normalization constant is simple to calculate. 
In this case, let us write Uk,,,fe t ,h simply as Uk^m^ fsee !4ap and we obtain 

lluk^m^H 2 = (A^.m./i) 2 "^ for m > 0, (19a) 
= (N khm , h ) 2 V, form = 0. (19b) 

So JV kll , 0> h = V- 1 ' 2 and JV k „, m> o, h = (2/V) 1 / 2 . 

For symmetric Lamb modes, from Eq. (|5b[) we calculate ki and k t as func- 
tions of fc||. The results are shown in Fig. [2] (a). As expected, for each value 
of fey, the wave- vectors ki and kt take only discrete values, but not as simple 
as the values corresponding to the h modes (|4aj). Each curve in Fig. [2] (a) cor- 
responds to a diffferent branch of the dispersion relation, u>ka,k t (,k»,m),s} where 
m = 0, 1, . . . denotes the branch number. Branches with bigger m are placed 
above branches with smaller m. 

As fcii increases, fc;(fcii,m) decreases and, after reaching the value 0, turns 
imaginary. On the other hand, fc t (fc||,m) first increases with fey, reaches a max- 
imal value and then decreases monotonically as fcii increases to infinity. Its 
decrease is bounded for all branches, except the lowest one, where, after reach- 
ing zero at some finite fcii , it turns imaginary in the lower-left quadrant of Fig. 
[2] (a). For the clarity of the calculations, the imaginary values of ki and kt are 
denoted as ini and in t , respectively, where ki and Kt take real, positive values. 
Both, K\ and Kt increase without limit as fcii increases to infinity. 
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Figure 2: The curves (ki(k\\),kt(k\\)) for symmetric (a) and antisymmetric (b) 
Lamb modes. On each curve k\\ varies from to oo, as we go from its right 
end to the left. The different curves correspond to different branches of the 
dispersion relation, i.e. to different values of fc; and k t at ku = 0. Both, ki and 
k t may take real and imaginary values. The imaginary values of ki, denoted by 
K\ are to the right of the vertical axis and the imaginary values of k t , denoted by 
K t are below the horizontal axis in both figures. As in the rest of the paper, d is 
the thickness of the plate. For these plots we used ci/c t = 1.66, corresponding 
to SiNa; plates. The + and — signs in the upper-right quadrants are only for 
orientation purposes. 



We encounter a similar situation for the antisymmetric Lamb modes (see Fig. 
[2]b). The only marked differences between the symmetric and antisymmetric 
modes are the following: the asymptotic values of (d/2) ■ k t (k\\ — > oo, m) are 
n, 2-7T, . . . for symmetric modes and 7r/2, 37r/2, . . . for antisymmetric modes; the 
maxima of (d/2) ■ fc t (fc||,m) are (2m + l)7r/2 for the symmetric modes and tyvk 
for the antisymmetric modes. Also, for the antisymmetric modes fc t (fc||,0) (i.e. 
the lowest branch) takes only imaginary values. 

Integrating |u s | 2 over the volume of the plate and using the transcendental 
equation (|5aj) . together with SnelPs law, u 2 — c 2 (k 2 + k 2 ) = cf(kf + k 2 ), we 
obtain the normalization constant for the symmetric modes in the quadrant I 
(upper-right) of Fig. [5] (a) 

(NlT 2 = A{4fc t 2 fc^cos 2 (M/2)((fcf + fc^)^-(fc ; 2 --^)^^) 

+ (fc t 2 - k 2 ) 2 cos 2 (M/2) {(k 2 t + k 2 /- + (k 2 t - k 2 )^M) 

+4k t kj(k 2 - kf) cos 2 (fc/d/2) sin(M) \ (20a) 
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in the quadrant II (upper-left), 
(iVf)- 2 = A|4fc t 2 fcfeos 2 (M/2)((^ + ^)^^-(^-fcf)0 

+ (fc t 2 - fc 2 ) 2 cosh 2 (^/2) ((fc 2 + fc 2 )^ + (fc 2 - * 2 )^g^) 

+4fc t fc 2 (fc t 2 - fc 2 ) cosh 2 {nid/2) sin(M) | (20b) 
and in the quadrant III (lower-left), 

(iV/'T 2 = A|4 K 2 fc 2 eosh 2 (^/2)^ 2 + fc 2 )^^-(. 2 -fc 2 )0 

+ ( K 2 + fc 2 ) 2 cosh 2 ( K( d/2) (% 2 + klfj^M + _ fc|)f) 
-4/t t fc 2 (/« 2 + fc 2 )cosh 2 (Kid/2)sinh(K t d) I (20c) 



Similarly, the normalization constants for symmetric modes in the three 
quadrants are 

(iVi)- 2 = A|4fc 2 fc 2 sin 2 (M/2)((fcf + ^)^ + (fc 2 -fcf)^^) 

+ (fc t 2 - fc 2 ) 2 sin 2 M/2) ((fc 2 + k\) d - - (fc 2 - fcf )^^) 



-4fc t fc 2 (fc t 2 - fc 2 ) sin 2 (fc ( d/2) sin(M) j (21a) 
(iVf)- 2 = A|4fc 2 fc 2 sin 2 M/2)^ 2 + fc 2 )^^ + ( K 2 -fc 2 )0 

+ (fc t 2 - fc 2 ) 2 sinh 2 (M/2) ((fc t 2 + fc 2 )^ - (fc 2 - *jf)^j^) 
-4fc t fc 2 (fc ( 2 - fc 2 ) sinh 2 (M/2) sin(M) 1 (21b) 



and 

(K 11 )- 2 = A{4. 2 fc 2 sinh 2 Kd/2)(( K 2 + fc 2 )^^ + ( K 2 -fc 2 ) 
+(«» + fc 2 ) 2 sinh 2 (M/2) ((. 2 + fcjj)!^M _ _ 
-4K 4 fcjj(K 2 + fc 2 )sinh 2 (K;d/2)sinh(K t d) I (21c) 



As one can notice, in the quadrants // and /// the trigonometric functions 
are replaced by hyperbolic functions because of the change fc; — > in quadrant 
//-and fc; — > iK;, fe< — ► iK t ~in quadrant 



Although Eqs. (f2"0")) and (f2"Tj) . with all their cases, are convenient because 
they are ready to use in practical calculations, we shall pack each set in one 
equation with the help of the complex wave-vector components k t = k t + int 
and ki = ki + iKi'. 

= A {aiS,! 3 ^ | coBC^d^)! 8 (d*,l a + *if)^^^ - (|*«| a - *f)^^^ 

+ n - *f n cos ( m/2)i 2 ( ( n 2 + + w _ 

-4fc 2 |cos(M/2)| 2 («t(N 2 + fef)sinh(Ktd) - h(\k t \ 2 - k 2 )sin(k t d) 



(22a) 

iV- = A {4|^| 2 A:{f| si n (A:^/2)| 2 + ^^^^h^) + (| ^ |2 _ ^^Ei^O" 

- k\ | 2 | sin(M/2)| 2 ((|^| 2 + fc 2 )5i^M - { \k t? - kjf-^fl 
m(k l d/2)\ 2 [n t (\k t \ 2 + k 2 ) sinh(/c t d) + k t (\k t \ 2 - fcjj)sin(M) 



-4fc 2 | sin 



(22b) 



Equations (|2"2"|) are very general. For the phonon modes (i.e. for real k\\), k t and 
ki are either real or imaginary and the normalization constants (|20[) and (|21[) 
can be extracted from (|22[) by taking the limit of the redundant component of 
fct and fc; going to zero. The normalization constants (|19p . ([20]) . and (f2Tj) are 
chosen so that ||uk,, k t ,<r\\ = 1> f° r an y ^'^ll an< ^ ^ ne nex t section the 

wave- functions will be multiplied by still another constant, which will give the 
right dimensions to the phonon field. 



4 Quantization of the elastic field 

For the quantization of the elastic field we start from the classical Hamiltonian: 
U = J (^ + S ^ S A . (23) 

where Sij are the components of the strain field, which is the symmetric gradient 
of the displacement field, SV, = (Vsu)y = (diUj + djUi)/2. The canonical 
variables are the field, u, and the conjugate momentum, tv = pu, which satisfy 
the Hamilton equations 

u = (24a) 

* = — . (24b) 
ou 

Equation (|24b[l is nothing but the dynamic equation {T]). 
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In the second quantization, u and 7r become the field operators, u and 7ir, 
respectively. If we denote by bt k and 6 k|l u t <r> the creation and annihilation 
operators of a phonon with quantum numbers k|| , k t and polarization a (in 
the notation that we used before), then the real displacement and generalized 
momentum field operators, u(r) = u^(r) and 7r(r) = 7r'(r), are 



fi W = [ f k||,fc t , CT ( r )&k ll ,fe t , CT +fk || ,fe tlCT (r)&] £||ifcti 



and 



Tr(r) 



kii ,k t ,rr 



P [KM,Ar)b^M.^ + K ll .k t A T )~ b l ll ,k t ,a 

k|i,fe t ,cr 



(25a) 



kii ,k t ,a 



tkii.fc 



tA V ) b 



kii ,kt ,<t -^kii ,fe 



t ,a(r)6: 



kii ,k t ,a 



(25b) 



where fk|i,fc tl o-( r ) = Ci*kn,k t ,a{ r ) an d C is a real constant which we shall deter- 
mine from the commutation relations of the b operators. In Eqs. (|25p we do 
not take ki as a summation variable, since this is either zero (for h fields), or 
is determined by k|| and k t (via Eq. (|5"a|) or Eq. (|5b[) for the symmetric and 
antisymmetric case, respectively). 

From Eqs. ([23)1 we can extract the operators b and ¥ in terms of u and tt. 
In order to do this, let's first note from Q that 

[uk||,fc t ,h(r)]* = u_ k|| ,fe t ,/t(r), (26a) 

whereas 

K,fe t , 8 W]*=U-k| l ,i t , s (r) or [u k||)fct>s (r)]* = -u_ k|l , fct , s (r), (26b) 
depending weather k t is real or imaginary, respectively, and 

[uk,|,fc tl o(r)]* =u_ k|| ,fe t ,a(r) or [u k|| jfct . Q (r)]* = -u_ k|| , fct , a (r), (26c) 
depending weather ki is real or imaginary, respectively. So 

/ f^A.aW-fk^fe^'W^r = a aM ,k l C 2 5 (T , IT .5^ t _^8 ktiK (27) 
J v 

where by f T we denote the transpose of the vector f and a ai k t ,ki = ac ~ 
cording to Eqs. (|26b|) and iggcj) , Multiplying ([2"5ajl and (|25b|) by f^t k _(r) and 
integrating over V, we get 



kti ,fc t ,iT 



(r)u(r) d 3 r 



fk lhfctlCT (r)^(r)d 3 r = 



-iW|,ti,a 



,c 2 



&k||,fe t ,cr a CT,fe t ,fc i fcLk |h fc t , ( j 



(28a) 
(28b) 
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Solving the system we obtain 



1 



I ^. kt Jr)u(r)d 3 r+— J— f ft (r)*(r)d 3 r 



(29) 



and obviously, 



St 



1 



ku.fct.o- 2C 2 



C ifct ^(r)u(r)d 3 r 



(30) 

Using the canonical commutation relations, [u(r),u(r')] = [7r(r), 7r(r')] = and 
[u(r), 7r(r')] = iK5(r— r'), we obtain the commutation relations for the operators 
b and w: 



and 



provided that 



[&k||,fc t ,<7,& k / fe / )ff ,] - (W^kj^fc^; 



c = 



2pW k ||,fc t , < r 



(31) 



Using Eqs. (|25[) . with the proper normalization of f , we can write U (|23p in 
operator form, 



/ ' /! '- X E [fk^^W^^^+fk^fc^WSk,,,^,. 

kii ,fc t ,cr k' ,k' er' 



f k ; t , k ' t ,k> r *> (r)b^ , Ktk ^ + f£, )fci|fc , (r)fot , >fc , jfc , ;CT , 
l - J d 3 r J2 [^[fk lh fc t , CT (i-)]^ k|l ,fc t , ( x + ^[fk||,fc t , CT (r)]^k l |,fc t . [ x 

kii ,k±,a k' ik'a' 



XC ijk l ■ [9fc[fk||,fc i ,<r(r)]l^ || ,fc t)(7 +9fcK || ,fe t)<T ('*)]/&k|| 1 fc t ,a 

E ^kn,*.,«r[^„ 1 fc t ,a^||.fct.«' + V 2 ] ■ 
kii ,fct ,f 



(32) 



As expected, the Hamiltonian of the elastic body can be written as a sum of 
Hamiltonians of harmonic oscillators. These oscillators are the phonon modes 
of the plate. 

We use this formalism elsewhere to describe the interaction of phonons with 
the disorder in amorphous materials [T31 [Tl] . 
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5 Conclusions 



The vibrational modes of a thin plate ((4]) are well known from elasticity theory 
[TT] . The purpose of the paper is to quantize the elastic field and for this we 
have to know if these modes, or part of them, form a complete set of orthogonal 
functions. But since the modes are the solutions of the eigenvalue-eigenvector 
problem of the operator L @, we showed that they form a complete set by 
proving that L is self-adjoint. 

Nevertheless, not all the functions of the form ((4]) are orthogonal to each- 
other, so to build the complete, orthogonal set of functions, we made use of 
a generic "momentum" operator, k|| = i(d x + d y ), which commutes with L. 
Since for a plate with infinite lateral extension or a finite rectangular plate with 
periodic boundary conditions at the edges the operator k|| is also self-adjoint, 
L and k|| admit a common, complete set of orthogonal eigenfunctions. The 
degenerate eigenvalues of k|| are, of course, the wave- vectors parallel to the plate 
surfaces, k|| , of real components. Therefore, the complete set of eigenfunctions 
are the ones given by Eqs. ([4|, with real k||. 

In Section [221 based on the hermiticity of the operator L, we showed that 
these functions ((4]) are indeed orthogonal to each-other and in Section [3] we 
calculated the normalization factors. 

Having all these ingredients, in Section [4] we presented the formal quanti- 
zation procedure which is applied elsewhere p~3|. [14] to calculate the thermal 
properties of ultra-thin plates at low temperatures, and to deduce some of the 
observed features of the standard tunneling model in bulk amorphous materials. 
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